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ABSTRACT 

We present deep 1.2 millimeter photometry of 37 stars in the young (5 Myr) Upper Scorpius OB 
association, sensitive to ~ 4 x 10~^ Mjup of cool millimeter dust. Disks around four low- and solar- 
mass stars are detected, as well as one debris disk around an intermediate mass star, with dust 
masses ranging from 3.6 x 10"'^ - 1.0 x 10~^Mj„p. The source with the most massive disk exhibits 
a transition-disk spectral energy distribution. Combining our results with previous studies, we find 
the millimeter-detection fraction of Class II sources has significantly decreased from younger ages, 
and comparison with near-infrared and Ha measurements indicates the present disks have undergone 
significant evolution in composition or structure at all radii. The disks of Upper Scorpius represent 
the tail-end of the depletion of primordial disks; while a few near-solar mass stars may still sustain 
giant planet formation, this process has finished around higher mass stars. 

Subject headings: circumstellar matter — open clusters and associations: individual (Upper Scorpius 
OBI) — planetary systems: protoplanetary disks — stars: pre-main-sequence 



1. INTRODUCTION 

Our understanding of circumstellar disks, the birth- 
place of planets, has been driven in the last decade largely 
by the Spitzer space telescope, which has allowed for 
the identification of young stars with warm, optically 
thick inner dust disks. Spitzer has observed all known 
nearby sites of recent star formation (e.g. c2d, FEPS, 
Gould's Belt Survey), revealing hundreds of new disk 
bearing systems with ages from 1 - 100 Myr. These 
disks have shown a remarkable variety of spectral energy 
distributions (SED), indicating the presence of large- 
excess primordial disks, weaker excess "anemic" disks, 
transition disks showing an emission gap at short wave- 
lengths, and debris disks showing weak emission charac- 
teristic of small amounts of second gener ation dust (for 
an overview, see IWilhams fc Cie^l20TTh . These many 
surveys have revealed a picture of circumstellar disk evo- 
lution in which inn er d isk dust dissipates on timescales 
of a few Myr (e.g. iHillenbrand .2006. ) . with inner disks 
around higher mass stars evolving more quickly (e.g. 
iCarpenter et all 120061 ). 

Optical and near-infrared spectroscopic surveys for ac- 
cretion indicators (e.g. Ha line equivalent width) have 
revealed a similar varie ty in accretion rates among young 
stars (e.g. 'Dahm|[2008D . The presence of accretion indi- 
cates the existence of gas in the inner disk. 

On the other hand, probing the outer disk, the main 
reservoir of dust, requires sensitive (sub-)millimeter pho- 
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tometry. At these wavelengths, the dust is optically 
thin, allowing for the direct measurement of dust mass. 
Past large mm-surveys in young (<10 Myr) regions have 
reached dust mass sensitivit ies of a few 10~^ M,7.,,p for 
standard assumptions (e .g. [Andrews fc Williams! 120051 : 
iRoccatagliata et aLll2009[ ). In this work, we present re- 
sults of a sensitive search for millimeter continuum in 
a sample of nearby stars in the Upper Scorpius region 
which have signs of significant evolution in the inner disk 
dust and gas. 

At a distance of 145 pc l|de Zeeuw et al.|[l999l ). Up- 
per Scorpius (Upper Sco) is one of the closest sites 
of recent star for mation. At an age of 5 Myr 
(jPreibisch fc Zinnec kcr 1999;), it has largely exhausted 
its interstellar gas and dust, reducing contamination 
from field emission. Spectroscopic s tudies indicate that 
few stars are stiU accre ting (e.g. iWalter et al.l 119941 : 
iDahm fc Carpenteill2009l ): this suggests that most disks 
have depleted their inner disk gas content. Spitzer stud- 
ies u sing IRAC, IRS, and MIPS (jCarpenter et al.|[2006l 
120091 ) directly revealed that the disks present in Up- 
per Sco are more evolved than disks in younger star 
forming regions. The majority of Upper Sco members 
(^80%), at all spectral types, exhibit no infrared ex- 
cess. SEDs for the infrared detections in Upper Sco 
indicate anemic or transition disks around K and M 
stars, and debris disks around stars across the entire 
sample (B through M spectral types). Similar results 

have been found in o ther 5 Myr associations (e .g. A 

Ori and NGC 236 2, 'Barr ado v Navascues et"al] 120071: 
IDahm fc Hillenbran d 2007). 

The variety of disk types around low mass stars in Up- 
per Sco suggest that an age of 5 Myr may be a time 
at which disk properties rapidly change, providing an 
ideal laboratory for examining the evolution of disks. 
Of all mechanisms affecting disk evolution, perhaps the 
one that has stirred the greatest interest in recent years 
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has been the possibihty of disruption by giant planets. 
Disk mass measurements wih provide direct constraints 
on the potential for on-going giant planet formation, in 
turn constraining planet formation models. 

We have carried out the most sensitive survey to date 
for millimeter emission from 5 Myr stars, reaching a dust 
mass sensitivity of ~4xl0~'^ Mj^p- In section [51 we de- 
scribe our target sample and our observations. In section 
[3l we discuss the inferred properties of our sources. In 
!|4l we discuss the properties of disks in Upper Sco, and 
in SjSl we compare to disks at younger ages and discuss 
the implications for planet formation. We summarize our 
findings in section [51 

2. SAMPLE, OBSERVATIONS, AND VALUES FROM THE 
LITERATURE 

We selected our sample from the 218 Upper Sco mem- 
bers o bserved in the near-infr ared (NIR) continuum sur- 
vey of [Carpenter et all (|2006[ ). They selected survey tar- 
gets that were identified as group members on the basis 
of stellar parameters so as to not bias towards disk bear- 
ing sources. Moreover, their targets were chosen so as to 
uniformly sample the stellar mass range in Upper Sco. 
We thus believe their sample to be the best representa- 
tion of the Upper Sco populatio n. Our sample focused 
on sources identified by Carpent er et al.l as having a NIR 
excess at 8 or 16 /im; we observed 8 of the 9 B and A 
(henceforth BA) excess sources, 20 out of 24 K and M 
stars (henceforth KM) with an 8 or 16 excess, and 3 
BA and 6 KM stars with no excess. 

Our targets were observed for continuum emission 
at 1.2 mm using the MAMB02 bolometer array 
(jKrevsa et al.l[l998l ) on the IRAM 30m telescope at Pico 
Veleta, Spain. Observations were conducted during the 
bolometer pools in 2008 (Nov. 13-18, 23) and 2009 
(February 24-25, March 8-9, and Nov. 8). Zenith opacity 
for our observations was typically ~0.2 - 0.3. Observa- 
tions were carried out to a target la sensitivity of 1 mJy, 
typically 20 minutes on source, in an ON-OFF pattern 
of 1 minute on target followed by 1 minute on sky, with a 
throw of 32". Observations during November 2009 were 
limited to 10 minutes on source, resulting in poorer sensi- 
tivity. In addition, deeper observations of objects show- 
ing 2(7 detections at the target RMS were carried out, 
resulting in the discovery of millimeter emission from a 
pair of sources (HIP 76310, ScoPMS 31) at a higher sen- 
sitivity than the median of the sample. Flux calibration 
was carried out using Mars, Saturn, or Jupiter (depend- 
ing on availability), and local pointing and secondary flux 
calibration was carried out using IRAS 16293-2422B. The 
data were reduced using the facility reduction software, 
MOPSICQ. 

We list the physical properties of the sources in this 
study in Table [TJ We include values for three objects 
from other studies: [PBB2002] USco J160823. 2-193001 
and [PBB20021 USco J160900.7-190852 (observed at 1.3 
mm in iCieza et ah , 20 08), and [PZ991 J161411.0 -230536 
(observed at 1.2 mm in lRoccatagliata et al.ll2009r). These 
are m embers of the parent sample from [Carpenter et al.l 
(|2006l ). and would have been observed as part of this 
work if photometry did not already exist in the literature. 

^ http: //www.iram.es /IRAMES / mainWiki/CookbookMopsic 



2.1. Spectral type and stellar mass 

The spectral types of BA stars are from 
iHernandez et al.l (|2005f ) , an d for KM st ar s we 
use spectral types from iPreibisch et al.l ()1998l 
20021). For BA stars , we adopt stellar masses from 



Hernandez et al.l (I2005D. who used the stellar isochrones 
of iPalla fc Stahleri (l2000l). For K M stars, we use masses 



from Kraus & Hillenbrand (2007*), who us ed the models 
of D 'Antona & Mazzitelh (1997) and iBaraffe et all 
(jl998l) for solar-mass and low mass stars, respectively. 
Two stars lack a mass estimate in the literature, 
[PBB2002] USco J160823.2-193001 and [PBB2002] USco 
J160900. 0-190836; we assigned these stars masses ac- 
cording to their spectral type, as in . Kraus & Hillcnbraii^ 
(|2007[ ). Comparison of the KM star m asses to other 
estimates in the literature (e.g. IPreibisch fc Zinnecked 
Il999f ) suggests a typical stellar mass uncertainty of 
10-20%. Both spectral type and estimated stellar mass 
are listed in Table [TJ 



2.2. Accretion 

The Ha equivalent wid th, W> (Hq;), was taken from 
iHernandcz et al. (2005^|Preib5ch ct al. (1998, 2002), 
iDahm fc Carpenteri (120091 ). and lRiaz et al.l (|2OO60 . This 
is then used to determine the accretion state of the 
stars. Traditionally, KM stars with WA(Ha) of lOA or 
greater were considered to be accreting, and referred 
to as Classical T-Tauri stars (CTTS). Those with 
lower or undetected emission were considered to be 
non-accreting, and referred to as Weak-line T-Tauri 
stars (WTTS). Recent work has refined this classi- 
fication system (iBarrado v Navascues fc Martini 120031 : 
iWhite fc Basril |2003| ). identifving spectral tvpe depen- 
dent W\{}la) that unambiguously indicate accretion. 
Here, we adopt the system of Whi te & Basril — WA(Ha) 
in emission of 3, 10, 20, and 40A indicate accretion for 
K0-K5, K7-M2.5, M3-M5.5, and M6 and later spectral 
types, respectively. While there are a few accreting sys- 
tems among the KM stars of Upper Sco, none of the BA 
stars in Upper Sco show Ha in emission, suggesting that 
none of these stars still experience accretion. 

We note that the use of an accretion indicator to 
demonstrate the presence of gas in the disk is subject 
to two main problems: accretion and photospheric Ha 
emission can be difficult to distinguish at low equiva- 
lent widths, and accretion may va ry by large amounts 
on timescales as short as days (jBouvier et al.l 120071 : 
iNguven et aIll2009D . Among the 20 KM stars m our sam- 
ple with multiple epochs of WA(Ha) measurements, two 
show values which would produce different classifications 
as CTTS or WTTS over a decade timescale. In these 
two cases ([PBB2002] J160545.4-202308 and [PBB2002] 
J160702. 1-201938), we classify the sources as CTTS. We 
list the WA(Ha) and accretion state of our sources in 
Table dl 

Our observations incl ude 8 of the 10 CTTS in the orig- 
inal sample of Carpente r et al.l (|2006l ). and 17 of the 95 
WTTS. There are' an additional 22 KM stars with no 
Ha measurement in the literature, of which we include 
1. The millimeter-detected sources we include from the 
literature are split between accreting and non-accreting 
systems; 1 is a CTTS, and 2 are WTTS. 
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2.3. Infrared disk classification 

Table [T] also shows the disk classification based on 
SEDs. SEDs were con structed using 2MASS J, H, and 
K band photometry (iSkrutskie et all I2006D. 4 .5. 8.0. 
and 16 /Ltm photometry from iCarpenter "etall (l2006h 
and 24 and 70 /Ltm photometry from ICarpenter et aLl 
([2009). We used extinction estimates froni the liter a- 
turc (Her nandez et all 120051 : iPreibisch craT 1998, 2002^ 
and the R=3.1 reddening law of ,Fitzp_atrick (1999) to 
deredden NIR photometry. 

For the B A stars in Upper Sco, w e use the SED clas- 
sifications of ICarpenter et all (I2009D. Among B A stars, 
they found only debris disks. ICarpenter et al.l classified 
stars in Upper Sco as hosting debris or primordial disks 
on the basis of the strength of the NIR excess and com- 
parison of NIR colors with those of disks in the younger 
IC348 region. 

For KM stars, the classification is based upon the slope 
of th e SED at NIR wavelengths, as in ([Greene et al.l 
119941). Assuming \Fx oc A", sources exhibiting a near 
infrared power law index greater than 0.3 are considered 
Class I, —0.3 < n < 0.3 are flat-spectrum, —1.6 < n < 
—0.3 are Class II, and n < —1.6 are Class III. The classes 
roughly trace the expected evolution of inner disk dust. 
We use the broadest baseline applicable to our entire 
sample, calculating the index from K-band and 24 fim 
photometry. 

We classify the K2 star [PZ99] J160421. 7-213028 as 
hosting a transition disk, based on its lack of NIR excess 
at wavelengths less than 16 iim and steeply rising excess 
at longer wavelengths (e.g. Cicza et al. 2007). 

The Upper Sco sample of ICarpenter et al.l (2006) con- 
tains 19 Class II sources, comprising 15% of the KM 
stars, and 107 Class III sources. We observed 16 Class II 
and 9 Class III sources. The three additional detections 
included from the literature are Class II sources. 

Due to the variety of disk types among KM stars, we 
summarize their disk classifications in Table [21 where we 
show the number of each type o f object in the parent 
sample of ICarpenter et all (|2006f ) (applying our Class II 
/ III and CTTS / WTTS criteria as described above), 
the number included in this work, and the number of 
millimeter detections. 

3. RESULTS AND ANALYSIS 

We detect 1 out of 11 BA stars, and 4 out of 26 
KM stars. Our median 3(t sensitivity is 2.8 mJy We 
report observation times, fluxes, and 3(T upper limits 
in Table [31 We also include fluxes from the literature 
for [PBB2002] USco J160823. 2-193001, [PBB2002] USco 
J160900.7-190852, and [PZ99] J161411.0-230536. 

Overall, millimeter observations exist for 9 of the 10 
CTTS in the original sample, of which 4 are detected. 
18 of the 95 WTTS have been observed for millimeter 
emission, with 3 detections, and the single millimeter- 
observed M star lacking Ha information is not detected. 

Examining disks by NIR classiflcation, 6 Class II 
sources are detected in millimeter emission. The addi- 
tional two millimeter-detections are the debris disk of 
HIP 76310, and the transition disk of [PZ99] J160421.7- 
213028. 

We show the SEDs of the 8 millimeter detected Up- 
per Sco sources in Figure [H We also display the upper 



quartile, median, and lower quartile SEDs o f disk bearing 
sourc es in the younger Taurus association ([Furlan et al.l 
[20061) for comparison. Most sources lie in the lower quar- 
tile of disk emission in Taurus, implying extensive evo- 
lution or settling of the dust. For the object [PBB2002] 
USco J161420. 3-190648, the photospheric model appears 
to be a poor fit at short wavelengths. This source is likely 
being viewed close to edge on. 

3.1. Mass estimates 

Dust masses were estimated from the 1.2 mm flux 
density assuming the dust to be optically thin, and are 
listed in Table [31 Under this assumption, the dust mass 
( Mdust) will b e direc tly proportional to the flux (F^), as 
in IHildebrandl ([1981 : 



Assuming a dust mass opacity {K,y) of 1 cm^/g at 1000 
GHz and opacity power law index /3=1 ([Beckwith et al.l 
IT990i) . as weU as a characteristic dust temperature (Tc) 
of 20 K (the median te mperature of Taurus disks; 
[Andrews fc Williamsl [2001 . we can then estimate dust 
masses as Mdust = 1-5 x 10"^fi.2mm(mJy) Mj^^at the 
145 pc distance of Upper Sco. For the two objects ob- 
served at 1.3 mm, [PBB2002] USco J160823. 2-193001 
and [PBB2002] USco J160900.7-190852, we must adopt 
a different constant: Mdust = 1.9 x 10^'^Fi.3mm(mJy) 
Mjup- While there is considerable theoretical uncer- 
tainty in these terms, the values are consistent with pre- 
vious disk studies. The addition of far-infrared photome- 
try will allow for detailed modeling of the dust structure 
of the disk and an improved dust mass estimate (e.g. 
iPinte et "aI1l2006f ). while multi-line sp ectroscopy is neces - 
sary to constrain the gas mass (e.g. iKamp et al.ll201l[ ). 
These are the subject of a future paper. 

Our median 3cr sensitivity of 2.8 mJy corresponds to 
an estimated dust mass sensitivity of 4.2x10"'^ Mjup- 
In order to estimate an upper limit to the average dust 
mass for non-detected BA and KM stars, we stack our 
non-detections. If many sources have mm-emission just 
below our 3(t limits, then we would expect stacking to 
result in a measurable flux. However, this procedure re- 
sults in a non-detection for both the stacked BA and KM 
observations. The 3cr upper limits of the stacked obser- 
vations, 0.89 mJy for BA stars, and 0.73 mJy for KM 
stars, do not differ from what would be expected from 
the noise in a single observation of the combined observa- 
tion times. These upper limits suggest that on average, 
BA and KM stars have dust masses less than 1.3x10"'^ 
and 1.1x10"'^ Mjup, respectively. 

4. DISCUSSION 

4.1. Comparison of detected and non-detected sources 

The millimeter-detected KM stars exhibit Ha emission 
similar to the non-detected KM stars. For each star, 
we calculate the mean Ha equivalent width from values 
in the literature, and show these in comparison to the 
dust masses in Figure [H We show sources with Ha in 
absorption as having an upper limit to emission of 0.1 
A. We calculate the Kendell's rank statistic, generalized 
for a censored data set ()Isobe et al..,1986) . for the KM 
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stars and find a Z-value of 0.028, indicating a 98% chance 
there is no correlation between the Ha equivalent widths 
of Upper Sco KM stars and their disk dust masses. Con- 
sidering just the mm-detected sources, there is ambigu- 
ous evidence for an anti-correlation (Z-value of 1.05, 29% 
chance of no correlation) . 

To carry out detailed comparison of the SED shapes, 
we plot the SED indices from 2.2 to 8 microns and 8 to 
24 microns in Fig. [3l defining the SED index ni_2 = 

log(^^)/log(^) as in iFurlan et all (poM ). The cen- 
tral cluster of points correspond to Class II sources, while 
points along the left hand side show no excess at 8 /im, 
indicative of a void inner region. However, some of these 
sources exhibit 24 /im excesses indicative of some amount 
of cooler dust which may either represent transition or 
debris disks. The majority of millimeter-detected sources 
in Upper Sco exhibit NIR SEDs typical of other Upper 
Sco Class II sources. 

While the NIR and Ha properties of the majority of 
millimeter detected sources do not stand out among Up- 
per Sco source, disk dust mass does appear as if it could 
be correlated with stellar mass. In Figure SI we show the 
dust mass as a function of stellar mass. The highest dust 
masses appear clustered around stars of 0.7-1.2 Mq, sug- 
gesting a stellar mass dependence in the evolution of the 
mass of millimeter-sized grains in disks. 

We calculate the Kendall's rank correlation statistic 
and find that the comparison is ambiguous when consid- 
ering the entire sample. With a Z-valuc of 0.87, there is a 
39% chance of rejecting the null-hypothesis that the disk 
mass and stellar mass are unrelated. Noting that the BA 
stars only host debris disks, we examine the stellar mass 
- dust mass relation for KM stars. In this case, the rela- 
tionship is clear; with a Z-value of 3.17, there is only a 
0.2% chance of disk dust mass not being related to the 
stellar mass. In Upper Sco, K stars tend to have disks 
with higher dust mass than do M stars. This could re- 
flect numerous factors; these higher mass stars may have 
emerged from their natal clouds with higher mass disks 
or some mechanism may be active in this mass range to 
preserve disks. 

The observation that millimeter detected sources are 
not outliers in their NIR and Ha emission suggests that 
the process leading to low millimeter emission is not nec- 
essarily tied to the loss of disk gas nor small warm grains. 
However, the concentration of high dust masses around 
near-solar mass stars suggests a stellar mass dependent 

mechanis m in dust evoluti on. This ex pands upon the 

result s of [Carpenter et al.l (|2006D and [Carpenter et al.l 
(|2009f ). who found evidence that small dust grains de- 
plete more rapidly around stars of greater than solar 
mass than around stars of near-solar mass or less; our 
results suggest millimeter-emitting grains are retained 
longer around near-solar mass stars. 

4.2. Effects of detectability 

The NIR and Ha classification systems we use above 
serve as indicators for the presence of dust and gas com- 
ponents of disks. While in younger groups the two prop- 
erties are strongly correlated, in Upper Sco the major- 
ity of Class II systems show no evidence for accretion. 
The statistics of our observations and data collected from 
the literature indicate that we are seeing the tail-end of 



the primordial disk stage. All disks hosted by BA stars 
appear to be debris disks. Only 22 out of 127 (17%) 
KM stars show any evidence for remaining primordial 
disk material (either NIR excess. Ha emission, or mm- 
emission), and only 7 (6%) have disks detected in mil- 
limeter photometry. 

These remaining primordial disks appear to be in the 
last stages of dissipation — their SEDs lie below the Tau- 
rus median, indicating significant inner disk dust evolu- 
tion or flattening of the disk; they generally have low 
Ha equivalent widths, suggesting depleted gas reservoirs 
for sustained accretion; and their low mm fluxes indicate 
that the mass in millimeter-sized and smaller grains has, 
for the most part, decreased to the edge of detectability 
with current facilities. 

One must keep in mind, however, that each of these 
measures reaches a different mass sensitivity. Because 
of this, one should exercise caution in interpreting the 
statistics regarding Class II / Class HI sources and CTTS 
/ WTTS sources. 

Our millimeter observations reach a dust mass sensi- 
tivity of 4 X 10~^Mjup, while near-IR observations are 
sensi tive to dust masses of ~ 10~^ Mj^p (e.g. lCieza et al.l 
I2007f ). These correspond to disk mass sensitivities of 0.4 
Mjup and 10~ ^ Mjuv, assuming the ISM gas-to-dust ra- 
tio of 100 fe.g. lBeckwith et al.|[l990l ). 

To estimate an equivalent disk mass sensitivity for Ha 
observations, we use the empir ical relation of a c cretio n 
rate as a function of time of iHartmann et al.l ()1998D . 
Modeling the accretion rate (M) as a power law in time 
(M (X t~^, with stellar age t and rj « 1.5 — 2.8), one 
can estimate the remaining disk gas mass as Mgas ~ 
(0.5 - 2) X M(t) X t. The WA(Ha) cutoffs for classi- 
fying systems as accreting roughly correspond to accre- 
tion rates of lO''^" Mq/jt (e.g. Eiahm 2008). Thus, for 
our Upper Sco sources with t — 5 Myr, WA(Ha) ob- 
servations provide an estimated disk mass sensitivity of 
~ (0.3 — 1) X 1O~^M0. While this estimate spans a large 
range of values, and neglects many of the details of ac- 
cretion (e.g. accretion variability, observation geometry), 
it illustrates that Ha observations are likely sensitive to 
disks of similar mass as those detected by our millimeter 
photometry. Both WA(Ha) and millimeter photometry 
based estimates of the disk mass are at least 2 orders of 
magnitude less sensitive than near-IR observations. 

Therefore, one should expect the presence of numerous 
systems exhibiting infrared excesses, but lacking one or 
both of detectable Ha emission and detectable millimeter 
emission. This is indeed the case in Upper Sco, with 
only 4 sources exhibiting all three disk tracers, while 80% 
of accreting sources and all millimeter detected sources 
exhibit an excess at 24 /im and shorter wavelengths. 

4.3. Notes on selected sources 

HIP 76310: The sole mm-detected A star in our 
sample, HIP 76310, has the brightest infrared excess 
of the BA stars in Upper Sco. However, its infrared 
to stel lar luminosity ra tio (L/ij/L* = 2.1 x 10^"*, 
iDahm fc CarpenteiH2009l ) falls within the observational 
definition of a debris disk {Lju/L^ < 10~^). Our as- 
sumption of 20K for the average dust temperature of 
primordial disks is based on KM stars in Taurus and 
does not apply to the case of a debris disk around an 
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A star. The SED fit by iDahm fc Carpenteii implies a 
much higher dust temperature of 122K and thus a cor- 
respondingly lower dust mass (a factor of ~8) . Both our 
initial mass estimate (3.6 x 10~^ Mj^p) and revised dust 
mass of ~ 5 X 10"'* Mj^p are several orders of magni- 
tu de higher than the low er limit of 2.6 x 10~^M,7„p found 
bv IDahm fc Carpenteii highlighting the high NIR opti- 
cal depth achieved by small quantities of dust, and the 
importance of (sub-) millimeter observations to directly 
constrain the dust mass. This dust mass also fits the 
evolutionary trend of de bris disk dust masses as a fu nc- 
tion of time as shown in iWilliams k, Andrews! ()2006[) . 

[PZ99] J160421. 7-213028: The K2 star [PZ99] 
J160421. 7-213028 stands out not only for having a 
uniquely steep infrared SED in Upper Sco, but also for 
having the largest millimeter emission. The SED shape 
for this object is typical of a transition disk, where dust 
grains have been cleared from the inner disk but are 
still abundant at larger radii. The strong millimeter 
emission implies a large dust mass (0.10 Mjup). How- 
ever, this source is not clearly accreting; the Ha equiv- 
alent width of -0.57 A suggests that some process has 
greatly reduced or cut off accretion to the stellar surface. 
By current disk evolution models, the two most likely 
candidate processes are the onset of photoevaporation 
([Alex ander et al..,2006,) or the formation of a planet large 
enoug h to reduce stellar accretion below our detection 
limits (|Alexander fc Armitagel[2007| ). Detailed examina- 
tion of this disk requires millimeter imaging, which we 
defer to a later paper. 

Binaries: ScoPMS 31 and [PZ99] J161411.0-230536 
are detected in millimeter continuum and exhibit pri- 
mordial type infrared excesses, yet have binary compan- 
ions a t projected separations of 84 and 32 AU, respec - 
tively (jKohler et al.l[2000l[Metchev fc Hillenbrandll2009[ ). 
There is ambiguity between observing a circumstellar or 
circumbinary disk; both cases allow for the presence of 
cold grains at large radii. While the question of the effect 
of binarity on disk mass and lifetime is worth consider- 
ing, our numbers here are too small to provide a useful 
study. 

5. DISK EVOLUTION 

We can investigate the evolution of observable disk 
properties from 1 to 5 Myr by comparing the properties 
of KM sourc es in Upper Sco with those in Taurus (mean 
age 1 Myr, [Andrews fc Wilhamj 120051 AW0 5), and in 
IC348 (mean age ~2.5 Myr; iLeeeLaU [mil LWCU). 
These studies are among the largest millimeter surveys at 
younger ages, for populations that have served as bench- 
marks in many disk studies. The Upper Sco and Tau- 
rus samples are readily comparable, as our study has 
similar dust mass limits to AW05 (4 x 10~^Mj„p ver- 
sus 6 X 10~'^ Mjup). Comparisons with IC348 (dust mass 
limit of 2 X 10~^Mjup) will necessarily focus on the high 
mass end of the dust mass distributions. We examine 
trends in disk mass as a function of disk accretion and 
NIR excess, and the implications for planet formation. 

5.1. Accretion status 

At younger ages accreting systems dominate the pop- 
ulation of sources with millimeter continuum detections. 
In Taurus, 67 out of 76 (88±11%) millimeter detected 



sources with Ha measurements in AW05 are accreting, 
while in IC348, all 9 mm-detected sources are accreting 
(LWCll). There is no such clear distinction in Upper 
Sco where, of the 7 millimeter detections, 4 are accreting 
and 3 are not. 

Furthermore, as noted in Section 14.11 the Ha distri- 
butions of millimeter detected and non-detected sources 
in Upper Sco are similar. This is a marked contrast to 
the younger IC348, where millimeter-detected disks are 
biased towards larger Ha equivalent widths than non- 
detected disks. 

Examining the detection statistics in greater detail, we 
find that the disk detection rate among CTTS in Upper 
Sco is 44 ± 22%, significantly lower than that in Taurus. 
The detection rate for WTTS is similar. In Table (H we 
show contingency tables constructed for detections and 
non-detections of CTTS and WTTS in the two regions. 
Using the Fisher exact test, we find for a two-tailed distri- 
bution the mm-detection rate among CTTS has ^ 0.25% 
chance of being independent of region; the apparent drop 
in mm-detection rate for CTTS we see in Upper Sco is 
significant. For WTTS, the two-tailed Fisher's exact test 
indicates that the mm-detection rate for WTTS remains 
the same between Taurus and Upper Sco. While the 
mm-detection rate for WTTS has changed little from 1 
to 5 Myr, the detection rate for CTTS has decreased 
markedly. 

5.2. Dust mass distribution 

The millimeter-detection rate for Class II sources drops 
significantly from 1 to 5 Myr. In Taurus, 86% of Class 
II sources have millimeter detections (AW05), while in 
Upper Sco that fraction has dropped to 32±13% (6/19). 
Furthermore, the median disk mass decreases. Taurus 
Class II disks have a median dust mass of 3 x 10~^Mj„p, 
while in Upper Sco the median dust mass of detected 
disks is 1.3 x 10~^Mj„p. The true median is likely to be 
lower, due to the high censoring rate in our data. 

The Class HI sources in Upper Sco are unlikely to host 
additional massive disks. We did not detect any of the 
9 Class HI disks in our survey, and the inferred upper 
limit to the detection rate of <11% is consistent with 
the detection rate for Class HI sources in AW05, 7±4%. 

As with the accretion state, we construct contingency 
tables of the detection statistics of Taurus and Upper Sco 
Class II and Class HI sources. We show the number of 
detections of each type of object in Taurus and Upper 
Sco in Table H As with CTTS, we find that the Class II 
sources are detected at a lower rate in Upper Sco than in 
Taurus. The detection rate for Class HI sources, on the 
other hand, is not clearly distinguishable. 

Due to the variation in exposure times in our observa- 
tions, and hence sensitivities, we use the Kaplan-Meier 
product limit estimator fo r randomly censored data sets 
^Feigelson fc NelsonI I1985D to generate the cumulative 
dust mass distribution functions, F{>MDust), for Class II 
sources in these regions (Fig. [5]). The dust mass distribu- 
tion of Upper Sco Class II sources is clearly different from 
that of Taurus Class II sources, with the upper quartile 
of Upper Sco sources overlapping the lowest quartile of 
Taurus sources; 27±10% of Upper Sco sources have disks 
with a dust mass greater than or equal to 4 x lO"'^ Mj„j,, 
compared to 89 ± 4% of those in Taurus. We make use 
of the Gehan, Peto-Prentice, and log-rank statistics to 
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compare the highly censored dust mass distributions of 
Upper Sco and IC 348; the distributions are different, 
with mild statistical significance (there is only a 2.5 -7% 
chance that the distributions are drawn from the same 
parent function). 

5.3. Disk dissipation timescale 

Having established that the dust mass distribution 
has changed from 1 to 5 Myr, we examine the rate at 
which the highest dust masses change. Comparing sim- 
ilarly constructed mass distributions, LWCll find that 
the fraction of IC348 disks with Mdust greater than 
2 X 10~^Mj„p is approximately equal to the fraction of 
Taurus disks with Mdust greater than 40 x lQ~^Mjup, 
suggesting a systematic factor of 20 decrease in the dust 
mass of Class II sources from 1 Myr to 2.5 Myr. If this 
trend were to continue, one would not expect to find any 
disks with millimeter continuum emission in Upper Sco. 
However, our results indicate dust masses do not decline 
so greatly from 2.5 to 5 Myr; the fraction of disks with 
Mdust > 2 X 10~^Mjup is similar between IC348 and 
Upper Sco, and the highest masses for Class II disks in 
Upper Sco are lower than masses of IC348 disks by only 
a factor of < 2. 

These results indicate that for Class II sources, the 
observed dust mass rapidly declines from 1 to 2.5 Myr, 
but may only slowly change from then until 5 Myr. In 
contrast, the fraction of disks exhibiting IR-excess con- 
sistent with primordial dust continues its marked de- 
cli ne. In Taurus, 50% of KM stars are Class I or Class 
II (jLuhman et al.ll2010l ). We include Class I sources here 
as their disks consist of primordial material. In IC348, 
only 39% of KM stars are Class II, w ith 1% of sources 
being Class I (Spit zer photometry of iLada et al.l 120061: 
iMuench et al.ll2007f ). In Upper Sco the Class II fraction 
is 15%, with no Class I sources. 

While there is no clea r evidence for a large age spread 
among Upper Sco stars. iPreibisch et al.l ()2002D could not 
rule out a 1-2 Myr age spread. This allows for the possi- 
bility of some overlap in stellar ages with the stars of IC 
348 (where there may be an age spread of several Myr). 
It is possible that the mm-detected sources in Upper Sco 
are also the youngest sources. If these sources were only 
marginally older than the sources observed in IC 348, the 
small decrease in the maximum mm-emission for Class 
II disks would be consistent with a continuation of the 
steep drop in mm-emission seen in comparing Taurus to 
IC 348. 

However, this argument could be equally applied to 
IC 348 and Taurus, with the corresponding claim that 
only the youngest disks are detected in the survey of 
LWCll, and the most massive disks in Taurus are also 
the youngest. Assuming similar age spreads in these re- 
gions, this scenario would not change the general picture 
we find of rapid change in the greatest disk masses over 
several Myr, followed by a more gradual decline. 

There have been recent attempts to trace the evolu- 
tion of disk properties as a functioii of st ellar age (e.g. 
Ilsella et al. I^OOOt IGuilloteau et al. |[20Tll) . The size of 
star forming regions, however, suggest distance uncer- 
tainties to individual stars large enough to cause 30% or 
greater uncertainties in luminosity estimates; other fac- 
tors (e.g. photometric errors) can make this uncertainty 



even higher. In turn, 30% luminosity uncertainties are 
large enough to cause a ge uncertainties of several Myr 
(|Hillenbrand et al.l [20081) . Until the age of individual 
stars can be accurately determined, interpretation may 
need to be restricted to consideration of average ages. 

Comparison of HR diagrams for these regions indicates 
that the median age in U pper Sco is older than T aurus 
and IC 348 (see Figure 1 of iHillenbrand et al.ll2008[ ). The 
comparison between Taurus and IC 348 is more ambigu- 
ous; across all temperatures, the median sequences could 
be equivalent within the luminosity errors. However, at 
temperatures of logT^ff less than 3.625 (corresponding 
to a spectral type of about K6), the median Taurus se- 
quence has luminosity equal to or higher than the IC 348 
sequence, suggesting that the median Taurus age for KM 
stars is likely younger than IC 348. These interpretations 
of the HR diagrams are consistent with the decrease in 
Class II and CTTS fraction seen from Taurus to IC 348, 
and then further to Upper Sco. While the uncertainty in 
absolute ages of star forming regions, age spread within 
those regions, and ages of individual stars complicate in- 
terpretation, the overall trend is clear. 

Significant progress has been made in recent years 
in modeling the dust evolution in disks. The amount 
of millimeter sized grains will decrease over time as 
dust aggregates to form centimeter and larger bod- 
ies, to which millimeter p ho tometry is iii s ensitive 
(jPullemond fc DominikI |2005| ). IBirnstiel et all (|2010bf ) 
modeled mm-slopes and fluxes in disks where millimeter- 
size grains remained present due to an equilibrium state 
between coagulation and fragmentation. Their models 
are consistent with obs ervations of young di sks at 0.5- 
2 Myr, and later work (jBirnstiel et al.ll201lD shows the 
same mechanisms can maintain populations of micron- 
sized grains that produce NIR excess emission. Observa- 
tions of millimeter continuum emission from older disks, 
such as this work, provide important constraints on the 
later time evolution of dust, and perhaps constrain the 
mechanisms leading to loss of observable dust (e.g. radial 
drift, agglomeration into centimeter and larger grains). 

In summary, there appears to be a rapid evolution of 
disks from 1 to 2.5 Myr, an interval during which both 
millimeter-dust masses and disk fractions rapidly drop. 
While the infrared disk fraction continues to decrease 
from 2.5 to 5 Myr, the dust masses of the remaining 
disks do not rapidly change. In addition, disks appear 
to reach a low baseline of dust mass where there is no 
longer a correlation between dust and inner disk gas. 

5.4. Implications for planet formation 

In addition to providing a point of comparison for the 
evolution of circumstellar disk mass as a function of age 
and inner disk properties, our disk dust mass measure- 
ments provide new insights on giant planet formation. 

Under the assump tion of a primordial 100:1 gas-to- 
dust mass ratio (e.g. iBeckwith et al.l [T990 [) as if planet 
formation has yet to begin. Upper Sco retains only a sin- 
gle known Minimum Mass Solar Nebula disk (MMSN, 10 
Mjup, WcidcnschiUing 1977), the 10 Mj„p disk around 
[PZ99] J160421. 7-213028. If planet formation were to 
begin now in Upper Sco, then the region could not form 
Jupiter mas s planets at t he frequency observed in the 
field (~10%.lJohn sonl2009f) . However, the disks in Upper 
Sco are clearly evolved, and the dust has grown substan- 
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tially. Jupiter mass planets must have already formed, or 
be well on their way to formation, as seen in the younger 
sources of IC348 (LWCll). The discovery of planetary 
companions to the Upper Sco me mbers IRXS J160929.1- 
210524 and GSC 06214-00210 (jLafreniere et all [25ol 
llreland et alllMlll serve as confirmation to this scenario. 
Due to our greater sensitivity, however, we can push this 
argument further to examine the formation of Neptune 
mass planets. 

At the 5 Myr age of Upper Sco, large rock and ice 
protoplanets may have formed, or may be in the pro- 
cess of formation. There are still uncertainties in the 
growth of grains from centimeter to planetesimal (10^ 
m) sizes, but simulations indicate growth from micron to 
centimeter size grains occurs in as little as 2 x 10^ years 
|Birnsticl et al. 2010a), and planetesimals can aggregate 
to form protoplanets with enough mass to accrete ga s 
(~0.03 Mjup) in less than 1 Myr (jPoUack et al.lll996D . 
Under core accretion models, the formation of a gas gi- 
ant planet requires that gas be retained in the disk until 
the accreting protoplanet reaches the critical mass (~0.1 
Mjup) at which runaway accreti on begins; this process 
could take from less than 1 Myr fAlibert et al."200_^ to 
several tens of Myr ([Pollack et a l. 1996), depending on 
several factors such as disk surface density, formation ra- 
dius, and the effects of planetary migration. 

Debris disks are expected to be seen in systems in the 
final stages of giant planet formation, as well as after 
the conclusion of giant planet formation. Our results 
verify that BA stars in Upper Sco have no residual cold 
millimeter-size dust above debris levels; these disks are 
unlikely to retain substantial gas mass, signaling the end 
of giant planet formation. 

However, some lower mass KM stars have NIR or mm 
continuum emission and / or Ha emission indicating the 
presence of gas and dust. The median inferred mass for 
the handful of mm-detected disks, under the assumption 
of the primordial 100:1 gas-to-dust ratio, is ~0.1 MMSN. 
By linear scaling of the MMSN, such a disk could pro- 
duce a planet approximately twice the mass of Neptune 
(Mjv = 0.05 Mjup) if planet formation were to begin 
now. Only 4 sources have disks with inferred masses of 
0.1 MMSN or greater, implying ^2% of stars in the 5 
Myr Upper Sco group retain disks able to form a 2Mn 
planet. In comparison, results of the first 4.5 months 
of the Kepler mission suggest that up to 19% of field 
stars have Neptune-s ize planet candidates within 0.5 AU 
([Borucki et al.l [2"011[ ). The formation of Neptune mass 
planets, in addition to the formation of Jupiter mass 
planets, must be at an advanced stage by 5 Myr. 

6. SUMMARY 

We have carried out the largest survey to date for mil- 
limeter continuum emission in a group of 5 Myr stars, 
an intermediate age in the evolution from ubiquitous gas 
and dust rich primordial disks at 1 Myr to, largely, gas 
and dust poor debris disks by 10 Myr. The overall picture 
of Upper Sco is of a region where a handful of primordial 
disks remain, with small disk masses, weak NIR excesses, 
and low Ha emission. We have found: 

1. 5 new 1.2 millimeter detections of 2.4 mJy and 
brighter, with estimated dust masses typically an 
order of magnitude lower than Taurus disks. 



2. The sole Minimum-Mass Solar Nebula equivalent 
mass disk is the transition disk of [PZ99] J160421.7- 
213028. 

3. The only BA star exhibiting mm-dust emission is 
the debris disk system with the highest 24 /xm ex- 
cess, HIP 76310. 

4. Other than [PZ99] J160421. 7-213028 and HIP 
76310, only Class II KM stars are detected in mm- 
emission. The Class II millimeter- fraction is 32%, 
a significant drop from 86% in Taurus. 

5. The mm-detection rate for CTTS is 44%, signif- 
icantly lower than the detection rate for CTTS 
in Taurus. Millimeter detections are roughly split 
among accreting and non- accreting sources, sug- 
gesting the action of mechanisms that can halt ac- 
cretion while preserving mm-sized grains. 

6. Dust mass in millimeter-emitting grains appears to 
decrease more slowly from 2.5 to 5 Myr than from 
1 to 2.5 Myr. 

7. We verify that BA stars with debris disk SEDs do 
not have large amounts of remnant cold dust, and 
find that with only four exceptions, KM stars with 
primordial disks lack sufficient mass to form new 
Neptune mass or larger planets. If the stars of Up- 
per Scorpius are to host giant planets at the same 
rate as field stars, then planet formation must be 
at an advanced stage. 
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Figure 1. SEDs of the millimeter detected sources discussed here, shown with Kurucz photospheric models l|Kuruc 3 [T979h and overlaid 
with the upper quartile, median, and lower quartile SEDs of Taurus KM stars (Furlan ct al. 2006) normalized to the H band flux. Red 
circles are dereddened fluxes from the literature (R=3.1. .Fitzpatrick |1999D . with open squares indicating original measurements. Our 
measurements are indicated as filled red squares. Inverted triangles indicate upper limits, la errors are smaller than symbols. 
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Figure 2. Estimated disk mass versus mean Ha equivalent widths for millimeter detected (red) and non-detected KM stars (gray) in 
Upper Sco. The displayed errors on M^^^^ only include flux errors; we neglect uncertainties due to adopted opacity and average disk 
temperature. There is no overall correlation between disk mass and W\{Ha). 
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Figure 3. Near- and mid-infrared SED indices, following the example of lFurlan et al.l 1)20061 ). Squares indicate BA sources, circles indicate 
KM sources, and plus signs indicate millimeter detected sources. The symbol color indicates accretion status, and "x" marks in the center 
and bottom left indicate the expected indices for a uniform flat disk and a stellar photosphere, respectively. The majority of detections 
have similar SED shapes to other Upper Sco Class II sources. 
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Figure 4. Estimated disk dust mass vs. stellar mass in Upper Sco (5 Myr). Solid circles indicate Class II or transition disks, solid squares 
indicate sources with debris disk-like IR excesses, and inverted triangles indicate non-excess sources. Millimeter detected sources are red 
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Figure 5. Cumulative dust mass distribution for Class II disks in Upper Sco, IC348, and Taurus, as well as Class III Taurus disks. The 
sample distributions are calculated using the Kaplan-Meier estimator to include 3a upper limits. The long vertical lines at the left indicate 
the lowest mass or upper limit in the sample. 
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Mathews et al. 



Table 1 

Target Properties 



Name 


RA 


Doc. 


Sp.Type^' 


WA(HQ)b 


SED 


Accretion 


M. 










[A] 


class 


state 





HIP 76310 
HIP 77815 
HIP 77911 
HIP 78099 
HIP 78996 
HIP 79156 
HIP 79410 
HIP 79439 
HIP 79878 
HIP 80088 
HIP 80130 
PBB2002 
PBB2002^ 
PBB2002' 
PBB2002' 
PBB2002^ 
PBB2002' 
PBB2002' 
PBB2002^ 
PBB2002' 
PBB2002' 
PBB2002' 
PBB2002^ 
PBB2002' 
PBB2002' 
PBB2002^ 
PBB2002' 
PBB2002' 
PBB2002^ 
PBB2002' 
PZ99 
PZ99 
PZ99 
PZ99 
PZ99 



USco 
USco 
USco 
USco 
USco 
USco 
USco 
USco 
USco 
USco 
USco 
USco 
USco 
USco 
USco 
USco 
USco 
USco 
USco 
J153557.8- 
J154413.4- 
J160108.0- 
J160357.6- 
J160421.7- 
RX J1600.7-2343 
ScoPMS 31 
[PBB2002J USco 
tPBB2002] USco 
[PZ99] J161411.0- 



J155624.8- 
J155706.4- 
J155729.9- 
J155829.8- 
J160210.9- 
J160245.4- 
J160357.9- 
J160525.5- 
J160532.1- 
J160545.4- 
J160600.6- 
J160622.8- 
J160702.1- 
J160827.5- 
J160900.0- 
J160953.6- 
J160959.4- 
J161115.3- 
J161420.3- 
-232405 
252258 
211318 
203105 
213028 



222555 
220606 
225843 
231007 
200749 
193037 
194210 
203539 
193315 
202308 
195711 
201124 
201938 
194904 
190836 
175446 
180009 
175721 
190648 



15:35:16.1 
15:53:21.9 
15:54:41.6 
15:56:47.9 
16:07:29.9 
16:09:20.9 
16:12:21.8 
16:12:44.1 
16:18:16.2 
16:20:50.2 
16:21:21.1 
15:56:24.8 
15:57:06.4 
15:57:29.9 
15:58:29.8 
16:02:11.0 
16:02:45.4 
16:03:57.9 
16:05:25.6 
16:05:32.2 
16:05:45.4 
16:06:00.6 
16:06:22.8 
16:07:02.1 
16:08:27.5 
16:09:00.0 
16:09:53.6 
16:09:59.3 
16:11:15.3 
16:14:20.3 
15:35:57.8 
15:44:13.3 
16:01:08.0 
16:03:57.7 
16:04:21.7 
16:00:44.6 
16:06:22.0 



-25:44:03.1 
-21:58:16.5 
-22:45:58.5 
-23:11:02.6 
-23:57:02.3 
-19:27:25.9 
-19:34:44.6 
-19:30:10.2 
-28:02:30.1 
-22:35:38.7 
-22:06:32.3 
-22:25:55.3 
-22:06:06.1 
-22:58:43.8 
-23:10:07.7 
-20:07:49.6 
-19:30:37.8 
-19:42:10.8 
-20:35:39.7 
-19:33:16.0 
-20:23:08.8 
-19:57:11.5 
-20:11:24.4 
-20:19:38.8 
-19:49:04.7 
-19:08:36.8 
-17:54:47.4 
-18:00:09.1 
-17:57:21.4 
-19:06:48.1 
-23:24:04.6 
-25:22:59.1 
-21:13:18.5 
-20:31:05.5 
-21:30:28.4 
-23:43:12.0 
-19:28:44.6 



J 160823.2- 
J160900.7- 
230536 



AOV 
A5V 
B9V 
AOV 
A9V 
AOV 
B9V 
B9V 
AOV 
A9V 
A9V 

M4 

M4 

M4 

M3 

M5 

M5 

M2 

M5 

M5 

M2 

M5 

M5 

M5 

M5 

M5 

M3 

M4 

Ml 

K5 

K3 

Ml 

MO 

K5 

K2 

M2 
M0.5V 
~K9 — 

K9 

KO 



< -0.1 



< - 

< - 

< - 

< - 

< - 

< - 

< - 

< - 

< - 
-5.40, 



0.1 
0.1 
0.1 
0.1 
0.1 
0.1 
0.1 
0.1 
0.1 
-5.51 



-3.60, -9.92 
-7.00, -6.91 
-250.00, -158.48 

-3.50 

-1.10 
-3.00, -2.70 



-6.10, 
-26.00, - 
-35.00, 
-7.50, 
-6.00, 
-30.00, 
-12.30, 
-15.40, 
-22.00, 
-4.00, 
-2.40, 
-52.00, 



-8.61 

152.12 

-2.04 

-4.11 

-3.12 

-8.30 

-14.56 

-12.88 

-22.23 

-4.41 

-4.47 

-43.72 



<-0.1 
-3.15 
-2.40 
-11.57 
-0.57 

-21.06 



debris 
none 
debris 
none 
debris 
debris 
debris 
debris 
debris 
debris 
none 

II 

II 

II 

II 

III 

III 

II 

III 

III 

II 

II 

II 

II 

III 

II 

II 

II 

II 

II 

III 

III 

III 

II 

transition 
III 
II 

— n 

II 
II 



No 

No 

No 

No 

No 

No 

No 

No 

No 

No 
WTTS 
WTTS 
WTTS 
CTTS 
WTTS 
WTTS 
WTTS 
WTTS 
CTTS 
CTTS 
WTTS 
WTTS 
CTTS 
WTTS 
WTTS 
CTTS 
WTTS 
WTTS 
CTTS 
WTTS 
WTTS 
WTTS 
CTTS 
WTTS 

CTTS 
WTTS 
CTTS 
WTTS 



2.20 
2.00 
2.80 
2.20 
1.80 
2.20 
2.50 
2.50 
2.30 
1.70 
1.90 
0.24 
0.24 
0.24 
0.36 
0.13 
0.13 
0.49 
0.13 
0.13 
0.49 
0.13 
0.13 
0.13 
0.13 
0.13 
0.36 
0.24 
0.60 
0.87 
0.99 
0.60 
0.68 
0.87 
1.12 
0.49 
0.64 



193001 
190852 



16:08:23.2 
16:09:00.8 
16:14:11.1 



-19:30:00.9 
-19:08:52.6 
-23:05:36.2 



-6.00, -2.75 
-12.70, -20.08 
0.96, 0.80, 0.38 



0.71 
0.71 
1.35 



" Spectral types from lHernandez et all l2005f l: [Preibisch et "aP 11199811200^ 

^ WA(Ha) fromlHernandez et all J2005I) : |P7eibisch et all Ill99a [200l): IDahm fc Carpenteil II2009I): IRiaz et aTl 120061) 



Table 2 

KM classification summary 



Disk class 




CTTS 




WTTS 


Unknown Ha 




COO 


This work 


mm 


COO 


This work 


mm 


C06 


This work mm 


Class II 


8 


8 


4 


11 


11 


2 






Transition 








1 


1 


1 






Class III 


2 


1 





83 


7 





22 


1 



Note. — Summary of KM star classifications, showing the number of each type in the parent 
sample jCarpenter et ani200^ ), the number in our study, and the number detected in millimeter 
continuum. 



Millimeter Survey of Upper Sco 



Table 3 

Time on source, measured flux, and 3cr Limits 









Name 


^souTce 

[minutes] 


F^(1.2mm) 
[m.Tvl 


M Oust 

\10~^ Ml. .J 


HIP 76310 


40 


2.41±0.61 


3.6 


HIP 77815 


30 


<2.00 


<3.0 


HIP 77911 


40 


<2.33 


<3.5 


HIP 78099 


20 


<2.80 


<4.2 


HIP 78996 


22 


<2.81 


<4.2 


HIP 79156 


32 


<2.18 


<3.3 


HIP 79410 


28 


<3.09 


<4.6 


HIP 79439 


34 


<2.63 


<3.9 


HIP 79878 


18 


<2.79 


<4.2 


HIP 80088 


28 


<2.47 


<3.7 


HIP 80130 


20 


<2.90 


<4.4 




PBB2002 


J155624.8-222555 


10 


<2.47 


<3.7 




PBB2002 


J155706.4-220606 


22 


<2.88 


<4.3 




PBB2002 


J155729.9-225843 


20 


<2.46 


<3.7 




PBB2002 


J155829.8-231007 


20 


<2.29 


<3.4 




PBB2002 


J160210.9-200749 


20 


<2.44 


<3.7 




PBB2002 


J160245.4-193037 


24 


<2.41 


<3.6 




PBB2002 


J160357.9-194210 


20 


<2.45 


<3.7 




PBB2002 


J160525.5-203539 


10 


<3.58 


<5.4 




PBB2002 


J160532.1-193315 


20 


<2.58 


<3.9 




PBB2002 


J160545.4-202308 


30 


5.15±0.76 


7.7 




PBB2002 


J160600.6-195711 


10 


<3.24 


<4.9 




PBB2002 


J160622.8-201124 


20 


<2.88 


<4.3 




PBB2002 


J160702.1-201938 


20 


<2.45 


<3.7 




PBB2002 


J160827.5-194904 


10 


<3.53 


<5.3 




PBB2002 


J160900.0-190836 


10 


<3.19 


<4.8 




PBB2002 


J160953.6-175446 


10 


<3.00 


<4.5 




PBB2002 


J160959.4-180009 


10 


<3.43 


<5.1 




PBB2002 


J161115.3-175721 


10 


<4.23 


<6 3 




PBB2002 


J161420.3-190648 


10 


9.79±1.14 


14.7 




PZ99 


J153557.8-232405 


20 


<2.78 


<4.2 




PZ99 


J154413. 4-252258 


20 


<2.78 


<4.2 




PZ99 


J160108. 0-211318 


20 


<2.68 


<4.0 




PZ99 


J160357.6-203105 


20 


<2.47 


<3.7 




PZ99 


J160421. 7-213028 


8 


67.51±1.44 


101.3 


RX J1600.7-2343 


20 


<2.50 


<3.8 


ScoPMS 31 


50 


2.70±0.61 


4.1 




PBB2002J USco J160823.2-193001 




21.90±4.70=' 


40.9 




PBB2002] USco J160900.7-190852 




12.30±2.00=' 


23.0 




PZ99] J161411. 0-230536 




3.5±l.l'= 


6.0 



^ 1.3 mm flux from lCieza et all 1 120081) . 

^ 1.2 mm flux froni fRoccatagliata et all j200g|) . 



Table 4 

Comparisons with Taurus disks 



1 CTTS 


WTTS 


Class II 


Class III 


mm-detected? 


yes 


no 


yes 


no 


yes 


no 


yes 


no 


Taurus 


67 


7 


9 


52 


64 


10 


4 


50 


Upper Sco 


4 


5 


3 


16 


6 


13 





9 


P(/Tau = fus) 


0.0025 


f» 1 


<<0.001 


f» 1 



Note. — Contingency tables for comparing detection statistics of CTTS / WTTS 
and Class II / Class III sources in Upper Sco and Taurus. P(/Tau — fus) shows the 
result of the two-tailed Fisher exact test for each table, showing the probability for 
the mm-detection rate in Taurus being equal to the min-detection rate in Upper Sco 
for each type of object. 



